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Abstract 

Control systems over networks with a finite data rate can be conveniently modeled as hybrid (impul- 
sive) systems. For the class of nonlinear systems in feedfoward form, we design a hybrid controller which 
guarantees stability, in spite of the measurement noise due to the quantization, and of an arbitrarily large 
delay which affects the communication channel. The rate at which feedback packets are transmitted from 
the sensors to the actuators is shown to be arbitrarily close to the infimal one. 

1 Introduction 

The problem of controlling systems under communication constraints has attracted much interest in recent 
years. In particular, many papers have investigated how to cope with the finite bandwidth of the commu- 
nication channel in the feedback loop. For the case of linear systems (cf. [2l |4l [30l [3 [Sj [23l [311 [3] to cite a 
few) the problem has been very well understood, and an elegant characterization of the minimal data rate - 
that is the minimal rate at which the measured information must be transmitted to the actuators - above 
which stabilization is always possible is available. Loosely speaking, the result shows that the minimal data 
rate is proportional to the inverse of the product of the unstable eigenvalues of the dynamic matrix of the 
system. Controlling using the minimal data rate is interesting not only from a theoretical point of view, but 
also from a practical one, even in the presence of communication channels with a large bandwidth. Indeed, 
having control techniques which employ a small number of bits to encode the feedback information implies 
for instance that the number of different tasks which can be simultaneously carried out is maximized, results 
in explicit procedures to convert the analog information provided by the sensors into the digital form which 
can be transmitted, and improves the performance of the system We refer the reader to |2S] for an 

excellent survey on the topic of control under data rate constraints. 

The problem for nonlinear systems has been investigated as well (cf. [Tll[Tni[211[lSl[ni[171[lH])- In [Ti]. 

the author extends the results of [5] on quantized control to nonlinear systems which are input-to- state 
stabilizable. For the same class, the paper [T^ shows that the approach in [3T] can be employed also for 
continuous-time nonlinear systems, although in |16| no attention is paid on the minimal data rate needed 
to achieve the result. In fact, if the requirement on the data rate is not strict, as it is implicitly assumed 
in [T6j, it is shown in [29] that the results of [16] actually hold for the much broader class of stabilizable 
systems. This observation is useful also to address the case in which only output feedback is transmitted 
through the channel. This is investigated in [28| for the class of uniformly completely observable systems. 
The paper |24| shows, among the other results, that a minimal data rate theorem for local stabilizability of 
nonlinear systems can be proven by focusing on linearized system. To the best of our knowledge, nan local 
results for the problem of minimal data rate stabilization of nonlinear systems are basically missing. Never- 
theless, the paper [27] has pointed out that, if one restricts the attention to the class of nonlinear feedforward 
systems, then it is possible to find the infimal data rate above which stabilizability is possible. We recall 
that feedforward systems represent a very important class of nonlinear systems, which has received much 
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attention in recent years (see e.g. [331 Uni [HI HHj j to cite a few), in which many physical systems fall 
and for which it is possible to design stabilizing control laws in spite of saturation on the actuators. When 
no communication channel is present in the feedback loop, a recent paper f|20j. see also |21| ) has shown 
that any feedforward nonlinear system can be stabilized regardless of an arbitrarily large delay affecting the 
control action. 

In this contribution, exploiting the results of [5D] , we show that the minimal data rate theorem of [57] holds 
when an arbitrarily large delay affects the channel (in [27j , instantaneous delivery through the channel of the 
feedback packets was assumed). Note that the communication channel not only introduces a delay, but also 
a quantization error and an impulsive behavior, since the packets of bits containing the feedback information 
are sent only at discrete times. Hence, the methods of |20| . which are studied for continuous-time delay sys- 
tems, can not be directly used to deal with impulsive delay systems in the presence of measurement errors. 
On the other hand, the approach of [27] must be modified to take into account the presence of the delay. 
To do this, we use different proof techniques than in |27| . Namely, while in the former paper the arguments 
were mainly trajectory-based, here we employ a mix of trajectory- and Lyapunov-based techniques. As a 
consequence, our result requires an appropriate redesign not only of the parameters in the feedback law of 
[20] . but also of the encoder and the decoder of [27]. Moreover, differently from ^27j, the parameters of the 
feedback control law are now explicitly computed (see Theorem [1] below). See [15] for another approach to 
control problems in the presence of delays and quantization. 

In the next section, we present some preliminary notions useful to formulate the problem. The main con- 
tribution is stated in Section [3] Building on the coordinate transformations of [32l [20] , we introduce in 
Section |4] a form for the closed-loop system which is convenient for the analysis discussed in Section [5] In 
the conclusions, it is emphasized how the proposed solution is also robust with respect to packet drop-out. 
The rest of the section summarizes the notation adopted throughout the paper. 

Notation. Given an integer 1 < i < i', the vector (a^, . . . , a^) e M''^'+^ will be succinctly denoted by the 
corresponding uppercase letter with index i, i.e. Ai. For i = 1, we will equivalently use the symbol Ai or 
simply a. li denotes the ix i identity matrix. Oixj (respectively, lixj) denotes an i x j matrix whose entries 
are all (respectively, 1). When only one index is present, it is meant that the matrix is a (row or column) 
vector. 

If X is a vector, |a;| denotes the standard Euclidean norm, i.e. |x'| = x^x, while |a;|oo denotes the infinity 
norm maxi<i<„ \xi\. The vector [x'^ IJ^Y' "^iH be simply denoted as {x^y). Z-|_ (respectively, R_|_) is the set 
of nonnegative integers (real numbers). R" is the positive orthant of M". A matrix M is said to be Schur 
stable if all its eigenvalues are strictly inside the unit circle. The symbol ||M|| denotes the induced matrix 
2-norm. 

The symbol sgn(a;), with x a scalar variable, denotes the sign function which is equal to 1 if a; > 0, if a; = 0, 
and equal to —1 otherwise. If x is an n-dimensional vector, then sgn(.T) is an n-dimensional vector whose 
ith component is given by sgn(a;i). Moreover, diag(x) is an n x n diagonal matrix whose element (i,i) is Xi. 
Given a vector-valued function of time x{-) : M", the symbol | |.t(-)| |oo denotes the supremum norm 

||2;(-)l|oo = supjgjj^ |2;(i)l- Moreover. x{i'^) represents the right limit limt^f+ a:(i). In the paper, two time 
scales are used, one denoted by the variable t in which the delay is 6, and the other one denoted by r, in which 
the delay is r. Depending on the time scale, the following two norms are used: ||a;t|| = sup_26i<?<o + 

11x^11 =SUp_2r<^<o|2;('' + Cr)|- 

The saturation function [20] a : R ^ R is an odd function such that < (j'{s) < 1 for all s e R, a{s) ~ 1 
for all s > 21/20, and cr(s) = s for all < s < 19/20. Furthermore, (Ji{s) = eia{s/ei), with Si a positive real 
number. The functions Pi,qi : R"^*+^ R are those functions such that [3^ 

/ . (n-i)! , , {n-i)\ 

with pi{qi{ai, . . . ,a„), . . . ,g„(a„)) = a;, qi{pi{ai, . . . ,a„), . . . ,_p„(a„)) = a^. 
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2 Preliminaries 



Consider a nonlinear system in feedforward form [331 1211 UHl HH] , that is a system of the form 

/ X2{t) + hi{X2it)) \ 



±{t) = f{x{t),u{t)) 



(t) + /i„_i(x„(t)) 

u{t) J 



(2) 



where Xi{t) E K, is the vector of state variables Xi{t),Xi+i{t), . . . ,Xn{t), u{t) G M, each function hi is 

C^, and there exists a positive real number M > such that for all z = 1, 2, . . . , n — 1, if |Xi+i|oo < 1, then 



MX,+i)\<M\X. 



l+l\ 



(3) 



We additionally assume that a bound on the compact set of initial conditions is available, namely a vector 



is known for which 



\x.,{to)\ <h 



1,2, 



(4) 



We investigate the problem of stabilizing the system above, when the measurements of the state variables 
travel through a communication channel. There are several ways to model the effect of the channel. In the 
present setting, wc assume that there exists a sequence of strictly increasing transmission times {tfcj/cez+i 
satisfying 

Tm < tk+i -tk<TM , k eZ+ (5) 

for some positive and known constants Tjn,T]\j, at which a packet of N(ti;) bits, encoding the feedback 
information, is transmitted. The packet is received at the other end of the channel 9 units of time later, 
namely at the times 9k := tk + 9. In problems of control under communication constraints, it is interesting 
to characterize how often the sensed information is transmitted to the actuators. In this contribution, as a 
measure of the data rate employed by the communication scheme we adopt the average data rate [31] defined 
as 



Rav = hmsupj,_ 



E 

i=0 



tk — to 



(6) 



where X]j=o ^i^j) the total number of bits transmitted during the time interval [to, tk]. An encoder carries 
out the conversion of the state variables into packets of bits. At each time tk, the encoder first samples the 
state vector to obtain x{tk), and then determines a vector y{tk) of symbols which can be transmitted through 
the channel. We recall below the encoder which has been proposed in [57] , inspired by [31] [TB] , and then 
propose a modification to handle the presence of the delay. 



2.1 Encoder in the delay-free case 

The encoder in [22 is as follows: 

lit) = 0„ t^tk 

= m+9£i^{t),m,m) (7) 

i{t+) = M{t) t = tk 

y{t) - sgn(x(t) - t = tk, 

where ^, £ is the encoder state, y is the feedback information transmitted through the channel, A is a 
Schur stable matrix, and g£{x,£,,£) = 4~^diag [sgn (x — The rationale behind the encoder (O is easily 
explained (we refer the reader to [3TJ[T6l|27] for more details). During the time at which no new feedback 
information is encoded, that is for t tk, the encoder tracks the state of the system. This explains the 
first equation of ([7|), which is a copy of the system ([2]). The (positive) values which appear in the entries 
of £ are the lengths of the edges of the quantization region, which is a cuboid with center ^. This vector 
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docs not change during continuous flow. At t = tk, when a new feedback information must be encoded and 
transmitted, the encoder sphts the cuboid into 2" subregions of equal volume, and select the subregion where 
the state x lies. The center of this subregion is taken as the new state ^ of the encoder (see the third equation 
in ([7])), and communicates this choice of the new state to the decoder through the symbol y. Note that each 
component of y takes value in {0,±1}, therefore y can be transmitted as a packet of bits of finite length. 
In particular, if is on the left of Xi then +1 is transmitted, if it is on the right, then —1 is transmitted. 
Finally, the size of the quantization region is updated, as it is shown by the fourth equation in ([7])- 
The system ([7]) is an impulsive system f[Tl 1121 and its behavior is as follows. At t = given an initial 
condition £,{to) i C-ito) , the updates of the encoder state and y(t^) of the output are obtained. 

The former update serves as initial condition for the continuous-time dynamics, and the state S^{t),i{t) is 
computed over the interval [io,ti]- At the endpoint of the interval, a new update S,{ti) t is obtained 

and the procedure can be iterated an infinite number of times to compute the solution ^{t),£(t) for all t. 
At the other end of the channel lies a decoder, which receives the packets y{tk), and reconstructs the state 
of the system. The decoder is very similar to the encoder. In fact, we have: 



(8) 



m = fm)Mt)) 

i'it) = On t^tk 

v{t+) = hv{t) t = tk 

with gv{y^ v) = 4~^dia,g{y)iy. The control law is the well-known nested saturated feedback controller 

U{t) = a{tp{t)) = " A^'^" (^l^"'^Y'^n{t)j + (Tn-l (pn-1 (^k""^ ^^An- 1 W , 



where 



+ A,-l(i))...) , 

+ ... + (71 (pi 

the functions pi are defined in ([1]), and the parameters k, L and the saturation levels Si of (Ti(r) = ei<7{r / £i) 
are to be designed (see Theorem [Ij. Note that we assume the decoder and the actuator to be co-located, 
and hence this control law is feasible. If the encoder and the decoder agree to set their initial conditions to 
the same value, then it is not hard to see ([17j) that ^(i) = ilj{t) and (.{t) = v{t) for all t. Moreover, one 
additionally proves that ^(t) is an asymptotically correct estimate of x{t), and the latter converges to zero 



2.2 Encoders for delayed channels 

When a delay affects the channel, both the encoder and the decoder, as well as the parameters of the 
controller, must be modified. As in |20[|32j . we shall adopt a /mear change of coordinates in which the control 
system takes a special form convenient for the analysis. Differently from [27] . this change of coordinates 
plays a role also in the encoding/decoding procedure. Denoted by $ the nonsingular matrix which defines 
the change of coordinates, and which we specify in detail in Sectional the functions gg, g-p which appear in 
the encoder and, respectively, the decoder arc modified as 

gsix, C, £) = (4$)-idiag [sgn {^x - C))] £ , gv{y, y) = (4$)-Miag(2/)i/ . 

Now, observe that the decoder docs not know the first state sample throughout the interval [1^,1^ + 0], and 
hence it can not provide any feedback control action. The control is therefore set to zero. As the successive 
samples y{tk) are all received at times 9k = + 0, the decoder becomes aware of the value o{ ^ 6 units of 
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(11) 



time later. Hence, the best one can expect is to reconstruct the value of ^(i — 6) (see Lemma [T] below), and 
to this purpose the following decoder is proposed: 

^{t+) = m + 9v{y{t-e),v{t)) (10) 

u{t+) = Kvit) t^6k 
u{t) = a(^(t)) . 

Observe that, as x{t) — f{x{t), a{il;{t))), and bearing in mind that ip{t) reconstructs ^{t — 9), the state x{t) 
is driven by a control law which is a function of £^{t — 6). This justifies the appearance of %l]{t ~ 6), rather 
than ip{t), in the first equation of the decoder above. 

We also need to modify the encoder. Indeed, as mentioned in the case with no delay, for the encoder to 
work correctly, the control law and hence ip{t), must be available to the encoder. To reconstruct this 
quantity, the following equations are added to the encoder ([7]): 

uj{t) = f{uj{t),a{Lu{t-0))) t^Ok 
u;{t+) ^ uj{t)+g£{x{t~e),^{t^e),l{t^e)) t^Ok. 

The initial conditions of the encoder and decoder arc set as 

uj{t) = 0„ , i e [00 - e, 0o] , ^{h) = 0„ , ^(io) = 2M , 

V'(i) ^On,te[eo-e, 6o] v{eo) = 2M , 

and, finally, the vector y which is transmitted through the channel takes the expression 

y{t) = sgn{<^{x(t) ~ m)) ■ 
Overall, the equations which describe the encoder are: 



(12) 



The following can be easily proven: 

Lemma 1 In the above setting, we have: (i) Lj{t) ~ ip{t) for all t > tQ, (ii) ^(t — 0) = iplt) and v{t — 9) = iit) 
for all t > 6q. 

Proof. First observe that the solution of ([2]), (fTO|) . ([T2|) exists for all t > to. Further, the thesis is trivially 
true for t e [to, Oq], by ([Til)- Now, for aU t ^ 9k, uj{t) - ipit) satisfies 

±[u;{t)^m] = f{u;{t),a{u;{t-e)))-f{m:c^mt~0))) (13) 
By definition of y, after the reset the initial conditions become 

uj{e+) = ^(0o) + (4a>)-idiag[sgn(a>(x(to)-e(io)))]^(to) 

= i;{eo) + m-^toHeo) ^ He+) 

As for t e [9o,6o + 0] the right-hand side of is equal to zero, surely we have that uj{t) = ip{t) over 
such time interval provided that no state reset occurs. Even when a state reset occurs (this is the case for 
instance if the delay is larger than the sampling interval) we can draw the same conclusion. In fact, before 



u;{t) 


- f{u{t),a{Lu{t-e))) 




m 


= fm.c^i^m 


t^tk 


m 


= o„ 


t ^ tk 


u;{t+) 


= oj(t)+ge(x{t-e),i{t-t 


^),£{t~0)) t = 0k 


at-') 


= m+9£{^{t),i{t),i{t)) 


t = tk 


i{t+) 


= M{t) 


t = tk 


y{t) 


= sgn(<i>(x(0 - e(0)) 


t^tk 
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the reset takes place, say at time t, w„(f) = ipn{t) for what we have just observed. But then, exactly as in 
(|14p . immediately after the reset the two variables are equal, and they continue to be the same up to time 
Oq + 9 because the growth of their difference is identically zero and other possible resets produce no effect as 
before. Iterating these arguments shows (i). For the second part of the thesis, observe that, by definition, 
^(tg) = i(;{to), and that ^(i^) = ipiOk) implies ^{t+) = i^iO^) exactly as in HH). Moreover, at t ^ 9k, 
d[i{t -0)- i){t)]ldt = 0, for di{t - 0)/dt depends on a{uj{t - 9)), while d->p(t)/dt depends on a{^j{t - 9)), 
and we have already proven that iu{t) = ip{t) for all t >tQ. The second part of (ii) is trivially true. 

As anticipated, the encoder and decoder we introduced above are such that the internal state of the 
former is exactly reconstructed from the internal state of the latter. Hence, in the analysis to come it is 
enough to focus on the equations describing the process and the decoder only. 

3 Main result 

The problem we tackle in this paper is, given any value of the delay 9, find the control ([9]) and the matrices 
A, $ in and (fTO|) which guarantee the state of the entire closed-loop system to converge to the origin. As 
recalled in the previous section, at times tk, the measured state is sampled, packed into a sequence of N(tk) 
bits, and fed back to the controller. In other words, the information flows from the sensors to the actuator 
with an average rate Rav given by ([6]). In this setting, it is therefore meaningful to formulate the problem 
of stabilizing the system while transmitting the minimal amount of feedback information per unit of time, 
that is using the minimal average data rate. The problem can be formally cast as follows: 

Definition. System ^ is semi-giobally asymptotically and locally exponentially stabilizable using an 
average data rate arbitrarily close to the infimal one if, for any £ G K", 9 > 0, R > 0, a controller ([9]), an 
encoder (fT2|) . a decoder pO)) . and initial conditions dH, pT|) . exist such that for the closed-loop system with 
state X := {x, lu, ^, £, ip, v): 

(i) There exist a compact set C containing the origin, and T > 0, such that X{t) e C for all t >T; 

(ii) For all t>T, for some positive real numbers k, 6, \X{t)\ < k\\XT\ \ exp(— (5(t — T)); 
(in) Rav < R- 

Remark. In the definition above we are slightly abusing the terminology, since for the sake of simplicity no 
requirement on Lyapunov (simple) stability is added. Nevertheless, it is not difficult to prove that the origin 
is a stable equilibrium point, in the sense that, for each e > 0, there exists > such that, ||a;eo|| < S^^ 
implies \x{t)\ < e for all t > 9q. 

Remark. Item (iii) explains what is meant by stabilizability using an average data rate arbitrarily close to 
the infimal one. As a matter of fact, (iii) requires that the average data rate can be made arbitrarily close 
to zero, which of course is the infimal data rate. It is "infimal" rather than "minimal" , because we could 
never stabilize an open- loop unstable system such as Q with a zero data rate (no feedback). 

Compared with the papers [33l [22l [TOl [18j , concerned with the stabilization problem of nonlinear feed- 
forward systems, the novelty here is due to the presence of impulses, quantization noise which affects the 
measurements and delays which affect the control action (on the other hand, we neglect parametric un- 
certainty, considered in |18j). In [34], it was shown robustness with respect to measurement errors for 
non-impulsive systems with no delay. In [20] , the input is delayed, but neither impulses nor measurement 
errors are present. Impulses and measurement errors are considered in [27], where the minimal data rate 
stabilization problem is solved, but instantaneous delivery of the packets is assumed. 
We state the main result of the paper: 

Theorem 1 // in (0) 

1 = 30£„ = 30^en-i = ... = 30"ei 

K > 6* max {6- {30)"+hi{n + l),8n^{8n{l + n'^)''-^ + 1)H , . 

\ K K Mk 1 ^^^^ 

< L < m"i|20-30"n4(n!)3' 8(1 n2)"-i V^n3(n!)3 ' (n + 1)!'^''^/ ' 

then system is semi- globally asymptotically and locally exponentially stabilizable with an average data 
rate arbitrarily close to the infimal one. 
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Remark. The parameters of the nested saturated controller are very similar to those in [20| . although k 
and L, e^'s arc, respectively, larger and smaller than the corresponding values of ^U\, to accommodate the 
presence of the quantization error. Moreover, in the proof, the values of the matrices A, $ are explicitly 
determined as well. 

Remark. This result can be viewed as a nonlinear generalization of the well-known data rate theorem for 
linear systems. Indeed, the linearization of the feedforward system at the origin is a chain of integrators, for 
which the minimal data rate theorem for linear continuous-time systems states that stabilizability is possible 
using an average data rate arbitrarily close to zero. 

4 Change of coordinates 

Building on the coordinate transformations in [20l [32] , we put the system composed of the process and the 
decoder in a special form (see ((23)) below). Before doing this, we recall that for feedforward systems encoders, 
decoders and controllers are designed in a recursive way [311 [331 [TOl HOI [57] . In particular, at each step 
i = 1,2, ... ,n, one focuses on the last n — i + 1 equations of system design the last n — i + 1 equations 
of the encoder and the decoder, the first i terms of the nested saturated controller, and then proceed to 
the next step, where the last n — i equations of ^ are considered. To this end, it is useful to introduce 
additional notation to denote these subsystems. In particular, for i ~ 1,2, ... ,n, we denote the last n — i + 1 
equations of ([2]) by 

/ x,+iit) + lM{X,+i{t)) \ 



V 



2;„(t) + /i„_i(X„(t)) 
uit) 



X,{t) = H,{X,+,{t),u{t)) = 
with u{t) = a{ip{t)), while for the last n ~ i + 1 equations of the decoder pI7)) we adopt the notation 



(16) 



*,(i) = H,i^,+iit),uit^9)) 

N,{t) = 0„_,+i 

= *,(t) + (4$O"'diag(y,(t-0))A^,(t) 

N,{t) = A,N,{t) 



(17) 



t = 



where Ni denotes the components from i to n of v. Moreover, for given positive constants L < M , k > 1, 
with M defined in ([3]), we define the non singular positive matrice^j as: 



1, • 



M 



M 



(18) 



where the functions pi are those defined in ([T]). Finally, let us also introduce the change of time scale 

t = Kr , (19) 

and the input coordinate change 

v{r) = Kpn(^^ti"-hi{Kr)^ . (20) 

Then we have: 

Lemma 2 Let i G {1, 2, . . . ,n} and 

t:^9/k, rk-.^tk/n, pk-=Ok/K.. (21) 



^The matrix $1 will be simply referred to as <&. 
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The change of coordinates il9\) . i20\) . and 



E,[r) = $,(*,(Kr)-X,(K(r-T))) 



(22) 



transforms system ll6\ )- [T7\ l into 





= r,z,(r) + 


ln^i+iv{r) 


+ ip,{Z,+i{r)) 






- T^E,{r)^ 


ip.,{E,+i{r) 
{r - r)) 


+ Z,+i(r-r)) 




Mr) 


= 0„_i+i 










- Z^{r) 








E,{r+) 


= E,{r)+A 


^^diag(sgn(- 


-E,{r)m{r) 




P^{r+) 


- A,P,(r) 






r ^ Pk , 



(23) 



whe 



1 1 
1 



1 1 
1 1 








1 




f^{Z^+l) 
fi+l{Zi+2) 

fn-l{Zn) 





whe 



provided that |Zi-|_i|oo < (AfK)/(L(n+ 1)!). 



P = n''^(7i!)^LK-i 



(24) 



(25) 



Proof. The proof is easily derived from analogous arguments in [201 [32]. However, for the convenience of the 
Readers, we are adding it in the Appendix. 

In the new coordinates ([T9|) - ([20|) . ([22]) . the controller ([9]) takes the form 

v[r) = -(T„(e„(r) + 2„(r - r) + cr„_i(e„_i(r) + z„_i(r - r)+ 
. . . + <7iiei{r) + Zi{r - t) + \i-i{r)) ...)), 

with Ai_i(r) = crj_i(ej_i(r) + 2j_i(r - r) + . . . + ai{ei{r) + zi(r - r)) . . .). 

5 Analysis 

In the previous sections, we have introduced the encoder, the decoder and the controller. In this section, in 
order to show the stability property, we carry out a step-by-step analysis, where at each step i, we consider 
the subsystem ([23]) in closed-loop with ([25]) . We first recall two lemmas which are at the basis of the iterative 
construction. The first one, which, in a different form, was basically given in [57], shows that the decoder 
asymptotically tracks the state of the process under a boundedness assumption. 

Lemma 3 Suppose is true. If for some i = 1, 2, . . . , n there exists a positive real number Z^+i such that 

E 

and, for all r > po, 
witl^ 



||^i+l(-)l|oo < 'Z'i+1 , 

lejWI <Pj{r)/2, j = 1 + 1,1 + 2,. ..,n, 
Pi+i {p~^) = Ai+iPi+i (p) p = Pk , 



■^The conditions are void for i = n. 

^In the statement, the continuous dynamics of the impulsive systems are trivial — the associated vector fields are identically 
zero - and hence omitted. 
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and Ai_|_i a Schur stable matrix, then for all r > po, 

|e^(r)| <K(r)/2, 
with pi(r~^) = pi[r)/2, for r — p^, if i = n, and 



1/2 * 



(26) 



r = Pk , 

if i & {1, 2, . . . , n — 1} , where * is a 1 x {n — i) row vector depending on Zi+i , £, and Tm- 

Proof. See |27| . For the convenience of the Readers, we arc adding the proof in the Appendix. 

The following observation is useful later on. As a consequence of the mean value theorem, it is not 
difhcult to realize (see also [U]) that, if ||2(-)||oo < Z, for some Z > 0, then e and p in (with « = 1) 
obey the equation^ 

e(r) = A{r)e{r) 

p{r) = 0,1 r ^ pk 

e[r^) = e(r) + 4^-^diag[sgn(— e(r))]p(r) 
p(r+) = Ap(r) r = Pk , 



(27) 



with 



A(r) 



ai2(r) ai3(r) 



023 (J^) 






ain-i[r) 
a2n-i(f) 






ainir) 
a2nir) 

an-inir) 




(28) 



and where the off-diagonal components of A, rather than as functions of {r, e{r), z(r — r)), are viewed as 
bounded (unknown) functions of r, whose absolute value can be assumed without loss of generality to be 
upper bounded by a positive constant depending on Z, £ and Tm- Now, concisely rewrite the system p7)) 
as([I]) 

e(r) = B{r)e{r) r ^ pk 



(29) 



with € = (e,p), |5fc(e)| > |e|/2. The following straightforward result shows that for the system above an 
exponential Lyapunov function exists (this fact was not pointed out in [27], where the proofs were not 
Lyapunov-based) : 



Corollary 1 There exists a function V{r,e) = V{r,e,p) : R+ x R" x 
and for all e = (e,p) G R" x R" for which |e| < |p|/2, satisfies 

ci|ep<V^(r,e)<C2|ep 

^ + ^B{r)e<~cM'' r + pk 
^(r+,gfc(e)) < l/(r,e) r = p^ 
dVir, e) 



such that, for all r G 



de 



< Ci\e\ 



for some positive constants ci, i = 1, 



Proof. By Lemma|31 for any r, |e(r)p + |p(r)p < 5|p(p^)p/4, where k is the index such that r S [pfc, Pk+i) 
Since the matrix A is Schur stable, there exist positive constants p, and A < 1 such that 

\p{pt)\ < pX''+'\pipo)\<pcxp{-\ir^po)Mpo)\ 



Again, we adopt the symbol A rather than Ai . 
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with A = I lnA|/Tm. Wc conclude that 

Hr)\ < ^/5/lpiexp{-\{r - po))\p{po)\ < exp{- X{r - po)) \e{po)\ , 

that is the solution is exponentially stable. Then, we can apply a standard converse Lyapunov theorem for 
impulsive systems, such as Theorem 15.2 in [T], and infer the thesis. 

The next statement, based on Lemma 10 in |20| . shows that a controller exists which guarantees the bounded- 
ness of the state variables, a property required in the latter result. For the proof, observe that the arguments 
of |20| hold even in the presence of a "measurement" disturbance e induced by the quantization, which can 
be possibly large during the transient but it is decaying to zero asymptotically. 



Lemma 4 Consider the system 

Z{r) = -ea 



-{Z(r-T)+e{r) + \{r)) 



where Z G M, e is a positive real number, and additionally: 



• A(-) and p,{-) are continuous functions for which positive real numbers A* and /i* exist such that, 
respectively, |A(r)| < A*, \tJ'{r)\ < H*, for all r > rp. 

• e(-) is a piecewise- continuous function for which a positive time r* and a positive number e* exist such 
that |e(r)| < e*, for all r > r^. 



If 



^^10' 12 



, A.g(o,^] ,e.e(o,^ 



< , and for all r > R, 



then there exist positive real numbers and R>0 such that \ \Z{- 

\Z{r)\ < 4(A» +e») . 

Proof See HO]- 

To illustrate the iterative analysis in a concise manner, the following is very useful (cf. the analogous 
inductive hypothesis in [20|): 



Inductive Hypothesis There exists > such that ||Zi(-)|| < Z^. Moreover, for each j 
\Gj{r)\ < pj{r)/2, for all r > po, and there exists Ri > t such that for all r > Ri, 



i,i + I, . 



1 



1 



\z,{r)\<-e,, |e,(r)|< — .— ^c, . 

Initial step (i — n) The initial step is trivially true, provided that r < 1/12, and e„ = 1/30. 
Inductive step The inductive step is summarized in the following result: 

Lemma 5 Let 

P<P,n< [20 • (30)"n]-i , T < T,„ < [6 • 30"+in(n + 1)]-^ . 
// the inductive hypothesis is true for some i G {2, . . . , n}, then it is also true for i — 1. 



(30) 



Proof. The proof only slightly differs from an analogous one in [20] in that wc have to deal with the presence 
of the quantization errors. The arguments arc as follows. If the induction hypothesis is satisfied for some 
i G {2, . . . , n}, then the control law v{r) in ((25|) becomes equal to 

-Wn{r) - ... - Wi{r) - CT,_i(w,_i(r) + K-2{'r)) = 

n 

-Xlt^j W -CT,_iK_i(r) + Aj_2(r)) , 

3=i 
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provided that \ej{r) + Zj{r — T)\ < ll£j/12 for each j = + . . . ,n, which is actually the case by hypothesis 
for all r > Ri + T. As a consequence, the components from i — 1 to n of the system in the variable z writes 
as 

Zi-i{r) = -<Tj_i(zj_i(r - t) + ei_i(r) + Xi-2{r)) + [z^(r) - Zi{r - r) + . . . + z„(r) 
-Zn{r - t)] - [e^(r) + . . . + e„(r)] + /,_i(Z,(r)) 
ze{r) = EjU+i^jW - EjU^jW - E"=»^j('^-'^)- (31) 
(^i-i{zi-i{r - t) + ei_i(r) + Xi-2{r)) + fe{Zi+i{r)) i=^i,...,n-l 

Znir) = - EjU ejir) ~ E"=. ^jX'^ " t") - cr^-i(zt-i(r - r) + e,_i(r) + Aj_2(0) 

As it is noted in |20j . from the expression of if(r), with £ — i,i+l,...,n, one proves that < 
{n + l)/4 + Pjnn/16, for < Pmu/IG. Hence, for each i ^ i,i + 1, . . . ,n, 

\zg{r) - Zi{r- t)| < I / ze{(T)da\ < r (—^— + 



so that 

'zi{r) - Zi{r - t) + . . . + Zn{r) - z„(r - t)] + /j_i(Zi(r))| < nr ( " ^ ''^ + Y^-f'™'^ 1 + ^^™'^ 



The following upper bounds easily follow from ([3^ 



1 ^ 1 n + 1 1 ^ 11 

— P,T,n < , Tn , rn ■ — Pmn < — • e,-i . 

16 320 4 ' 16 2 320 

By the inductive hypothesis, we additionally have |ei(r) + . . . + e„(r)| < ei_i/160. The bounds above show 
that in the ii_i(r) equation we have 

\[z,{r) - z,{r - r) + . . . + z„(r) - z„(r - r)] - [e,(r) + . . . + e„(r)] + /,_i(Z,(r))| < e,_i/80 . 

Furthermore, as a consequence of Lemma |3] and the inductive hypothesis, ei_i(r) satisfies |ei_i(r)| < 
Pi_i(r)/2, and in finite time becomes smaller than ei-ijiji ■ 160) < ei_i/30. Hence, it is possible to apply 
LemmaSl and infer that ||Zi_i(-)||oo < -^i-i and there exists a time Ri-\ such that |zi_i(r)| < ei_i/4 for 
ah r >Ri-i. 

Applying this lemma repeatedly, one concludes that, after a finite time, the state converges to the linear 
operation region for all the saturation functions, it remains therein for all the subsequent times, and hence 
the closed- loop system starts evolving according to the equations (recall ([21])) 



i(r) = rz(r) + Az(r — t) + Ae(r) + 'y9(z(r)) 
e(r) = B{r)e{r) r ^ pk 

z(?-+) = z{r) 

= 5fc(eW) r = pk, 

where: 

(i) r, A are matrices for which there exist q ~ {1 + n^)"^^, a ~ n, and Q = > such that 

(r + A)^g + Q(r + A) < -/, 

with IIQII < q and ||r||, ||A|| < a; 

(ii) There exists 7 > such that ip{z{r)) := [fi{Z2{r)) . . . /r,_i(Z„(r)) 0]-^ satisfies 

\'piz)\<l\z\ ; 

(iii) e = (e,p). 

In [501 the authors investigate the stability property of 

i(r) ~ Tz{r) + Az{r — t) + (p(z{r)) , 



(32) 
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that is the first component of system (j32p . with e = and no impulses. In the present case, e is due to the 
quantization noise and drives the z-subsystcm. The "driver" subsystem is described by the e equations of 
P2p . Hence, we have to study the stabihty of a cascade system with impulses. The Corollary [T] points out 
that there exists an exponential Lyapunov function for the e subsystem of (|32p . Based on this function, one 
can build a Lyapunov-Krasowskii functional to show that the origin is exponentially stable for the entire 
cascade impulsive system p2p . thus extending Lemma 11 in |20| in the following way: 

Lemma 6 Consider system iSS^) . for which the conditions (i)-(iii) hold. If 

1 , 1 

7 < — and r < 



8a2(8a(7+ 1)2 ' 

then, for all r > po, for some positive real numbers /i, 5, we have 

|(z(r),e(r))|<HI(^,e)pol|expM(r-po)). 
Proof. Following [20], we notice first that Vi{z) = z^Qz satisfies, by the requirement on 7, 

l>i < -i|z|2 + ^!i5^i±ll!r r \z{s)\'ds + 8q^a^T f \e{s)\'ds + 8q'a'\e\^ 

2 4 Jr-2T Jr-T 

The condition on t additionally gives 



2' ' 32 
The Lyapunov-Krasowskii functional 



Vi < -l\zf + ^r \z{s)\'ds+ I \e{s)fds + 8q'a'\e\' 

r-2r 



V2ir,z,e)^Viir,z) + ^ I I \z{l)\''dids + 2 I I \e{e)\'' d£ ds , 

r—2T J s J r — 2t J s 



16 

satisfies for r ^ pk 



V, < -i|zp + l/' \z{s)\^ds+[ |e(.)|2d5 + 8gV|ep 
-1 r W)\HI+\t\z\'^2 f |e(s)|2d£ + 4r|ep, 

-•-0 Jr-2T ° Jr~2T 



that is 

V2 < \z{s)\^ds^ I \eis)\'ds + 8{q'a' + lM 



32 



r-2r 



since t < 2. Now, recalling that |e(r)| < |p(r)|/2 for all r as a consequence of Lemma |3l and choosing the 
function V3{r, z, e) = V2(r, z, e) + XV{r, e) (see also [H] for a similar choice), with A > [1 + 8{q^a^ + l)]/c3, 
it is readily proven that 

V3 < -\\^\'-\^\'-^ f Ms)fds~ r Hs)fds, (33) 

4 6Z Jr-2T Ji — 2t 

for all r ^ p^. The inequality V3(r+, z, (/^(e)) < V^{r,z,e), for r = pk, also descends straightforwardly 
(cf. [Ej). Hence, we can conclude that the impulsive delay system is exponentially stable. Indeed, ai|(z,e)p < 
V3{r, z, e) < a2||(z, e)r|P, with ai = \m{Q) + Aci, 02 = Xm{Q) + 3x^/4 + Ac2. Moreover, 



V3 < -03 



AA/(Q)kr + Ac2|e|" + - / |z(s)|"rfs + 2T / |e(s)|"ds 

r-2T Jr-2T 



< -03^3 , 



with 03 ^ = max {4Aa/(Q), Ac2, 4t}. Over the interval [pk, Pk+i), we have ([6]) 

V'3(Pfc+i,z,fffc+i(e)) < V3{pk+i,z,e) < exp{-a3{pk+i - pk jjVsip^ , z, gk{e)) 
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A simple iterative arguments yields that for all r > po, 

ai|(z(r),e(r))p < V3(r,z,e) < exp{-a3{r ~ po))V3{po, z{po), e{pQ)) < a2||(z, e)pjp exp(-a3(r - po)) , 

that is 1(^(7'), e(''))| < y^a2/ai\\{z, e)p„\\ exp(— a3(r — pq)/2), which shows exponential stability of the impul- 
sive delay system. 

We can now state the following stability result for the system (|25p : 
Proposition 1 // the saturation levels Si are as in U5\). 

M K 

and 

0<r<r„i = [max|6 • (30)"+i?i(n + l),8n2(8n(l + + 

(35) 

0<P<P™ = [max{20-30"n,8(l + n2)"'iV"}]^^ 

then the closed-loop system h23\). with i = 1, and |ei(po)| £ PiiPo)/'^ for all i = 1,2, ... ,n, is semi-globally 
asymptotically and locally exponentially stabilized by the controller \25\) . 

Proof. Bearing in mind and that < 1. for i = 1,2, . . .,n, and (Mk) /{L{n + 1)!) > 1, then 7 in (ii) 
after ([32]) is equal to ^/nP, and the condition P < [8(1 + n^)"-^^"^ in actually implies 7 < l/{8q). 
Analogously, one can check that the remaining requirements on r and P in (|35p imply that all the conditions 
in Lemma [5] and [5] are true. These lemma allow us to infer the thesis. Indeed, by Lemma [5l the state after 
a finite time starts evolving in the linear operation region for all the saturation functions. Hence, from a 
certain time on, the state exponentially converges to zero, again by Lemma [S] 

Proof of Theorem Ql' The proof of the main result of the paper simply amounts to rephrase the proposition 
above in the original coordinates. Details are straightforward and we omit them. Only observe that the 
parameter k is chosen so as to satisfy the requirement on r in ([35|) . Once k is fixed, L is chosen to satisfy 
([M)) and the second inequality in ([55)1 . The equations of encoder and the decoder are introduced in and, 
respectively, ([TU)) . The matrices A and $ appearing in and, respectively, ([TU)) are designed in Lemma [3] 
and, respectively, defined in (fT8)) . 

As far as the minimality of the data rate is concerned, we proceed as in |27j . Observe that, by definition, 
Rav — '^n/Tm, and to guarantee Rav < R, we must have T„i > 2n/R. Since R can be arbitrarily close 
to zero, T,n can be arbitrarily large. Nevertheless, the stability result shown above holds for any value of 
Tm, for Tm (and hence Tm > Tm) affects the off-diagonal entries of A(r) and A only, and the exponential 
stability of the e equations (and therefore of system ((321) ') remains true, as it is evident from Lemma [3] (see 
also [27]). Hence, no matter what R is, we can always stabilize the system with Rav < R- 



6 Conclusion 

We have shown that minimal data rate stabilization of nonlinear systems is possible even when the com- 
munication channel is affected by an arbitrarily large transmission delay. The system has been modeled as 
the feedback interconnection of a couple of impulsive nonlinear control systems with the delay affecting the 
feedback loop. In suitable coordinates, the closed-loop system turns out to be described by a cascade of 
impulsive delay nonlinear control systems, and semi-global asymptotic plus local exponential stability has 
been shown. The proof relies, among other things, on the design of a Lyapunov-Krasowskii functional for 
an appropriate cascade impulsive time-delay system. If the encoder is endowed with a device able to detect 
abrupt changes in the rate of growth of Xn, or if a dedicated channel is available to inform the encoder about 
the transmission delays, then it is not difficult to derive the same kind of stability result for the case when 
the delays are time- varying and upper-bounded by 6. Similarly, by adjusting Tm in (O, it is possible to show 
that the solution proposed in this paper is also robust with respect to packet drop-outs, provided that the 
maximal number of consecutive drop-outs is available. The same kind of approach appears to be suitable 
for other problems of control over communication channel with finite data rate, delays and packet drop-out. 
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APPENDIX 

This appendix is added for the convenience of the Readers. 

Proof of Lemma [2] 

It is shown in gD] that, ^ and Zi{r) = ^,X,{nr) transforms ^ into 

Z,(r) = F,;(Z,+i(r),«(r)) 

EU+i^,{r)+v{r)+f.{Z..+i{r)) 



EU+2^Ar)+v{r)+f.+i{Z;+2{r)) 
v{r) 



(36) 
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Clearly, the equation is equal to the first equation of . Bearing in mind , and by differentiating 
Ei defined in (|22p . we obtain 



(37) 



E,{r) = F,iE,+i{r) + Z,+i{r-T),v{r-T))-F,iZ,+i{r-T),v{r-T)) 

= r,(£;,(r) + Z,(r - r)) + l„_.+iz;(r - t) + ^,{E,+i{r) + Z,+i(r - r)) 

-rj(Zj(r - r)) - l„„i+iw(r - t) - (y9j(Zj+i(r - t)) 
= TiEi{r) + ipi{E^+i{r) + Zi+i{r - t)) - (pi{Zi+i{r - t)) 

for r =/= pk, while for r ^ pk, we have: 

E,{p+) = <I>,(vI/,(«;p+)-X,(Ac(p,-T)+)) 

= M'^f^{9k) + {i<^^)-^dmgiY,{t+))N,{ek)-X,{tk)) (38) 
= $.(*,(0fe) - X,(<fe)) + 4-idiag(sgn($4X,(tfc) - E,{tk)]))m{9k) 
= $.(*,(0fe) - X.itk)) + 4-idiag(sgn($,[X,(tfc) - ^^{9k)]))N,{9k) , 

where the last equality descends from (ii) in Lemma [1] and implies 

E,{p+) = EM + 4-Miag(sgn(-S,(p,)))7V,(0fe) . (39) 

The thesis then follows if we observe that the variable Pi defined in satisfies 



P,(r+) = A,P,(r) r = pk. 
Proof of Lemma [3] 



(40) 



Recall first the equations in ((23|) which define the evolution of E^. Furthermore, by pT|) . the definition 
of $, and (HI, |ej(po)| < Pj{po)/2 for j = i, i + 1, . . . , n. For i = n, as |e„(po)| < Pn{po)/2, it is immediately 
seen that 

|e„(/9([)| = |e„(/3o) + 4~^sgn(-e„(po))Pn(po)| < 4"V«(po) 

which proves that |e„(/9(j")| < p„(pJ)/2, provided that A„ = 1/2. As e„(r) = 0, then |e„(r)| < p„(p(|)/2 for 
r G [po,Pi)- As Pn{r) = 0, also |e„(pi)| < p„(pi)/2, and iterative arguments prove that |e„(r)| < p„(p^)/2 
on each interval [pk, Pk+i)- Notice that the single trivial eigenvalue of A„ is strictly less than the unity. The 
equation for a in (|23|) writes as: 



( In-,: + 









ei{r) = ln-iEi+i{r) + ipi{Et+i{r) + Zi+i{r ~ t)) - ipi{Zi+i{7^ ~ t)) 

with a{r) E [0, 1] for all r. As both E'i+i and Zi+i are bounded, it is not hard to see [27] that there exists a 
positive real number Fi depending on Z^+i and i, such that, for r G [pk, Pfc+i), 

n 

e»(r) < e^ip-^) + F,(pfe+i - pfc) ^ Vjipt)/'^ , 
with |ei(/9(^)| < pi{pn)/4:. By iteration, the thesis is inferred provided that 

Pi{pt) = ip,(pfc) + p,rA/i„-,Ai+iPi+i(pfc) 



1 

> ::Pi{Pk) + F,{pk+i ~ Pk) Pj^Pk) 



]=i+i 



Note that, by the definition of Pi{Pk) above, Pi{Pk) ~ ^iPiiPk), with A^ the matrix in ([26]) . that shows A, 
to be a Schur stable matrix provided that so is A^+i. 
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